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ABSTRACT 

Context. A good correlation has been found in star-forming galaxies, between the soft X-ray and the far infrared or radio luminosities. 
The soft X-ray emission in star-forming regions is driven by the heating of the diffuse interstellar medium, and by the mechanical 
energy released by stellar winds and supernova explosions, both directly linked to the strength of the star formation episode. 
Aims. We analyze the relation between the soft X-ray and far infrared luminosities as predicted by evolutionary population synthesis 
models, aiming first to test the validity of the soft X-ray luminosity as a star formation rate estimator, using the already known 
calibration of the FIR luminosity as a proxy, and second to propose a calibration based on the predictions of evolutionary synthesis 
models. 

Methods. We have computed the soft X-ray and far infrared luminosities expected for a massive starburst as a function of evolutionary 
state, the efficiency of the conversion of mechanical energy into soft X-ray luminosity, the star formation history (instantaneous or 
extended) and dust abundance, and we have compared these predictions with observational values for 62 star-forming galaxies taken 
from the literature. 

Results. The observational Lsoftx/ipiR ratios are consistent with the model predictions under realistic assumptions (young starbursts, 
and efficiency in the re-processing of mechanical energy of a few percent), confirming the correlation between the diffuse soft X-ray 
emission and the star formation episode. 

Conclusions. The soft X-ray emission of the diffuse, extended gas surrounding massive star-forming regions, can be used as a star 
formation rate tracer. The empirical calibrations presented in the literature are supported by the predictions of evolutionary synthesis 
models, and by the analysis of a larger number of star-forming galaxies The calibrations are, however, biased towards galaxies 
dominated by relatively unevolved starbursts. 

Key words. ISM: kinematics and dynamics - ISM: supernova remnants - galaxies: starburst ~ X-rays: galaxies - X-rays: ISM 



.1. Introduction have proven to be invaluable for statistical studies of the star for- 

^ . mation history of the Universe. 
. . 1 The onset of massive star formation episodes in galaxies drives 

their observational properties in almost any wavelength range. Star-forming regions are also the source of conspicuous X- 

The UV and optical become dominated by the continuum of ""ay emission, generated by individual stars, by the injection 

massive, hot and young stars, as well as by the presence of neb- of l^S^ amounts of mechanical energy heatmg the mterstellar 

ular emission lines. After a few Myr of evolution, red super- medium, by supernova remnants, and by binary systems transfer- 

giant stars contribute to most of the neai- infrared emission. The ring mass to a compact primary d Cervino, Mas-Hesse, & Kunth| 

heating of interstellar dust particles by the powerful UV pho- M; [Persic & Rephaeb| | 2002|) . All of these individual compo- 

tons induces the thermal re-emission of large amounts of energy "^"1^ ^^e m prmciple directly hnked to the strength of the star 

in the mid and far infrared domain. The injection of ionizing formation episode, so that the X-ray luminosity could also be 

photons into the surrounding gas generates thermal radio emis- ^^^^ ^" estimator of star formation rates (SFR). 

sion, which is replaced by non-thermal emission as the ionizing Several authors have in recent years discussed the fea- 

power of the burst declines and the more massive stars begin to sibility of using the X-ray luminosity as an SFR estimator 

explode as supernovae. The direct relation between the strength Fabbiano & Shapley (2002) already concluded from the analysis 

of the star formation episode and the intensity of the different of 234 SO/a-IiT galaxies observed with Einstein, that the corre- 

observable parameters has allowed a number of star formation lation they found between the X-ray and the FIR luminosities in 

rate calibrators to be defined, such as UV continuum, emis- Sc-Iit ga laxies was due to the young stella r populations in these 

sion lines intensity, far infrared or radio luminosities (Ke nnicutti objects. Ranalli. Comastri. & Settil (l2003h proposed an empiri- 

ll998tlRosa-Gonzalez et a l. 2002; Bell 2003.) . These calibrators cal calibration of the soft (0.5-2.0 keV) and hard (2-10 keV) 

X-ray luminosities, based on their correlation with the far in- 

Send offprint requests to: J.M. Mas-Hesse frared (FIR) and radio luminosities, and using the known calibra- 



2 



Mas-Hesse et al.: Soft X-ray to Far Infrared luminosities ratio in star-forming galaxies 




15 20 25 30 " 5 10 15 2( 

Time (Myr) Time (Myr) 



Fig. 1. Evolution of LpiR (left) and Lsoftx (right) for IB (dashed line, right axis) and EB (solid line, left axis) models. IB models 
predictions are shown normalized to 1 Mq of gas transformed into stars, while for EB the luminosities are scaled to SFR - 1 Mq 
yr"'. LpiR has been plotted (from top to bottom) for E(B-V) = 1.0, 0.5 and 0.1. Lsoftx has been computed (from top to bottom) for 
exeif = 0.1, 0.05 and 0.01. 



tions of these parameters as proxies. [Grimm et al.l (l2003h stud- 
ied the correlation between the number of high-mass X-ray bi- 
naries (HMXB) and the SFR, deriving different calibrations of 
the hard X-ray lum inosity for low and high star formation rates. 
iPersic etaTI (|2004) obtained a calibration of the hard X-ray lu- 
minosity as a SFR estimator by assuming that most of the emis- 
sion in this range is associated with HMXB, an d scaling from 
the number of HMXB to the SFR of our Galaxy. iGilfanov et al] 
(|2004) confirmed the calibration of the hard X-ray luminosity as- 
sociated with HMXB, using slightly different slopes at high and 
low st ar formation rates. In a recent paper, iPersic & Rephaeli 
( |2007|) found indeed that the collective hard X-ray emission of 
young point sources correlates linearly with the star formation 
rate de rived from the far infrared luminosity. Strickland et al. 
(l2004ah demonstrated that the luminosity of diffuse X-ray emis- 
sion in star-forming galaxies is directly proportional to the rate 
of mechanical energy injection from the young, massive stars 
into the interstella r medium of the ho st galaxies. A similar re- 
sult was found by iGrimes et al.l (12005) from the analysis of a 
sample of starburst galaxies of different types (from dwarf star- 
bursts to ultraluminous infrared galaxies), which concluded that 
the mechanism producing the diffuse X-ray emission in the dif- 
ferent types of starbursts was powered by the mechanical energy 
injected by stellar winds and supernova e into the surrounding 
medium. Recently iRosa-Gonzalez et aTl (l2007l) confirmed the 
reliability of the soft X-ray luminosity as an SFR estimator from 
the analysis of a sample of star-forming galaxies in the Chandra 
Deep Field South at redshifts z - 0.01 - 0.67, using the UV 
continuum luminosity from GALEX as a proxy. 

While hard X-ray emission from star-forming regions may 
be dominated by binary systems, diffuse soft X-ray emission is 
generated by reprocessed, mechanical energy from stellar winds 
and supernovae explosions. This mechanical energy is related to 
the strength of the burst of star-formation, and can be calculated 
using evolutionary population synthesis models. In this paper, 
we analyze the correlation between the soft X-ray and FIR lumi- 
nosities in star-forming regions, both predicted by evolutionary 
synthesis models. We study the Lsoftx/^FiR ratio as a function 
of the evolutionary state, the efficiency of the conversion of me- 
chanical energy into soft X-ray luminosity, the star formation 
history (instantaneous or extended), and the dust abundance, and 
compare the computed values with observations taken from the 
literature. Our objective is to derive a calibration of Lsoftx as a 



tracer of the star formation rate, based on the predictions of evo- 
lutionary synthesis models, and to test the validity of the emp ir- 
ical calibration proposed bvl Ranalli. Comastri. & Settil (l2003h . 

In Sect. 2 we describe the evolutionary synthesis models that 
we use in the present study, in Sect. 3 we present the observa- 
tional data taken from the literature, and in Sect. 4 we discuss 
the predictions and the comparison with the observational val- 
ues. Throughout this work we have assumed Ho - 73 km s"' 
Mpc"' to convert fluxes into luminosities. 



2. Evolutionary synthesis models 

We have computed the predicted Lsoftx and LpiR val- 
ue s using the evolutionary p o pulatio n synthesis models 
of ICe rvifio. Mas-Hesse. & Kuiiflil (120021) (hereaf ter CMHK02 
model o, which are based on th e m odels of lArnault et al.l 
( 198%, iMas-Hesse & Kunflil (|1991|) and ICervino & Mas-Hessd 
(1994)- These models compute the evolution of a cluster of mas- 
sive stars, formed at the same time (Instantaneous Bursts, IB), 
or during an extended period of time (typically several tens of 
Myr) at a constant rate (Extended Bursts, EB), assuming differ- 
ent metallicities and Initial Mass Function (IMF) slopes. The 
CMHK02 models were developed to compute the evolution of a 
starburst during the first 30 Myr, following the onset of the star 
formation episode. Once the structure (number of stars of each 
spectral type and luminosity class at a given evolutionary time) 
of the stellar population is derived, the models are used to com- 
pute a number of observable parameters, from X-ray to radio 
wavelengths. 

In general, we assume a solar metallicity Zq, and a Salpeter 
IMF (0(m) ~ rrr'^-^^) with masses ranging between 2 Mq and 
120 Mq. The output of the models are normalized to the star 
formation rate (mass of gas transformed into stars per unit time, 
MQjr ') for EB cases, or star formation strength (mass of gas 
transformed into stars at the onset of the burst, in units of Mq) 
for IB scenarios. In both cases, the mass normalization corre- 
sponds to the mass integrated between 2 Mq and 120 Mq as- 
suming a Salpeter IMF. We emphasize that this normalization 
might differ significantly if other mass limits are considered. For 
example, the ratio between our mass normalization (2 - I2OM0) 
and the one assumed by .Kennicutt (1998) (0.1 - IOOMq) is 
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Fig. 2. Evolution of the Lsoftx/^FiR ratio computed for exeff - 0.01, 0.1 and E(B-V) = 0.1, 1.0. Left: predictions for IB models; right: 
predictions for EB models. 



M2^°/Mq°° = 0.293. It is important to remark that this is the nor- 
malization implicitly assumed when the SFR calibrations pro- 
posed by Kennicutt ( 1998) are used. Under our assumptions, the 
soft X-ray (0.4 - 2.4 keV) and far infrared luminosities were 
computed during the first 30 Myr after the onset of the massive 
starburst episode. 

One of the main sources of soft X-ray emission in a star- 
forming region is the diffuse gas heated by mechanical en- 
ergy from a starburst (massive stellar winds or supernova ex- 
plosions) into the surrounding medium. Its contribution is mod- 
eled by a Raymond-Smith thermal plasma with kT - 0.5 keV, 
which controls the fraction of the mechanical energy that heats 
the gas to X-ray temperatures, exeff. The models include the 
soft X-ray radiation emitted during the adiabatic phase of the 
Supernova Remnant (SNR), which is modeled by a compos- 
ite Raymond-Smith plasma with kT - 0.23, 0.76 and 1.29 
keV. The total energy emitted by the SNR, during the adia- 
batic phase, has been subtracted from the energy of each super- 
nova explosion when computing the injection of mechanical en- 
ergy. A more detailed d escription of our models is provided in 
ICervifio. Mas-Hesse. & Kunth (2002) . 

The contribution of the stellar atmospheres to the soft X- 
ray emission was neglected because it is expected to be two 
orders of magnitude lower than the emission from the diffuse 
gas. The contribution of high-mass X-ray binaries to the soft 
X-ray emission has in addition been neglected in this work. 
iMas-Hesse & Cerviflol d 19991) discussed the properties of the 
HMXB population expected to form during a massive star- 
formation burst. Binary systems become X-ray emitters when 
the primary collapses into a black hole or neutron star, the at- 
mosphere of the secondary has started to expand, and the sec- 
ondary is sufficiently close to the collapsed primary for mass 
transfer to begin. Mass is accreted onto the surface of the com- 
pact object and emits X-rays with a typical Lx ~ lO"^*^ erg s~', 
peaki ng at energies between 5 and 10 keV dPersic & RephaeUl 
I2OO2I) . The number of HMXB in a young starburst is depen- 
denton many free parameters. Following Mas-Hesse & Cerviflol 
(Il999l) . we estimate that only a few HMXB should be active after 
the first 5-6 Myr of evolution of starbursts that have transformed 
approximately 10* Moof gas into stars. This HMXB population 
should contribute a few times lO-'** erg s"' to the total X-ray lu- 
minosity, and a small fraction of this radiation to the soft X-ray 
emission. In all cases, the total contribution of HMXB, to the 
soft X-ray emission, remains 15% for IB and 10% for EB. Only 



if one or a few of these HMXB develop int o an ultraluminous X - 
ray source (ULX), with Lx ~ 10"° erg s ' dMiniutti et al.ll2006l) . 
could the X-ray emission, from soft to hard X-rays, of the entire 
galaxy, be dominated by the emission of HMXB, compared to 
that of the diffuse gas. 

Concerning the FIR emission, a thermal equilibrium of dust 
is assumed, implying that all the energy absorbed by dust, mostly 
originating from the UV continuum of the massive stars, is 
reemitted in the FIR range. In this paper, LpiR refers to the to- 
tal mid- and far-infrared luminosity integrated over the wave- 
length range 1 - 1000 |/m. We remark that this parametrization 
of LpiR implies a value of Lfir that is larger than that ca lculated 
using the FIR parameter propos ed bvlHelou et al.l (0*9881) . which 
is widely used in the literature. iHelou et al.l ( Il988l) showed that 
the FIR luminosity computed from the IRAS fluxes at 60 and 
100 yum would intercept about 70% (0.15 dex) of the total FIR 
luminosity from 1 to 1000 fim, assuming that there is a single, 
dominating component, at a temperature of 30 to 50 K. The dis- 
crepancy would be even larger in the presence of an additional 
warm dust component. This estimate of Lfir is c onsistent with 
the range (8 - 1000 //m) considered bv iKennicuttI ([l998) for the 
calibration of Lfir as an SFR estimator, since most of the FIR 
luminosity in starburst galaxies is emitted at wavelengths in the 
ra nge 10 - 120 jj.m. The models apply a Galactic extinction law 
(ICardelli et al.lll989 ) to the synthetic spectral energy distribu- 
tions, which is parametrized by the colour excess E(B-V). The 
energy absorbed due to extinction is calculated using the mod- 
els. It is assumed that all energy absorbed is reemitted thermally 
by the dust, within the mid to far infrared domain, i.e., within 
the ~ 8 - 1000 /im range. The models do not predict the shape of 
the FIR emission, since we do not make any assumption about 
the expected dust temperature. The presence of completely- 
obscured stars is not taken into account. Furthermore, the models 
assume that a fraction 1 - / of Lyman continuum photons is di- 
rectly ab sorbed by the dust, an d does not contribute to the ioniza- 
tion (see MezgeretaD ('1974). iMezged ([T978i) for the Galactic 
and iDegioia-Eastwood (1993) for the Large Magellanic Cloud 
HII regions, derived 1 - / values in the range 0.3 - 0.4. We have 
assumed (1 - /) - 0.3 i n this work, as proposed for starburst 
galaxies bv lBelfortet'an(ll987l) . 

The evolution of both Lfir and Ljoftx predicted by the mod- 
els are shown in Fig.lT] while their ratio Lsoax/Lfir is presented 
in Fig. |2] For IB models, the luminosities are shown scaled to 
1 Mq of gas transformed into stars at the onset of the starburst. 
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Fig. 3. Evolution of the Lsoftx/i^FiR ratio as a function of metal- 
licity and star formation regime. E(B-V)- 1 .0 and fxeff - 0.05. 



In the case of EB models, the luminosities are normalized to 
SFR = 1 MQjr LpiR is presented in the plot, computed for 
E(B-V) = 0.1, 0.5 and 1.0. It can be seen that FIR emission sat- 
urates rapidly for E(B-V) values above 0.5. In the remainder of 
this work, we consider the value of Lhr calculated by assum- 
ing that E(B-V) - 1.0, which can be considered an upper limit 
to the expected FIR luminosity. In the case of EB models, LpiR 
reaches an asymptotic value after approximately 10-15 Myr of 
evolution, when an equilibrium is reached between the number 
of massive stars that die, and those forming continuously. For 
coeval starbursts, LpiR declines rapidly after the first 5 Myr of 
evolution, when the most massive stars begin to end their life- 
times and stop heating the interstellar dust. 

isoftx is shown in Fig. [1] computed for fxeff values of 1%, 
5% and 10%. During the first few Myr, there is a rapid in- 
crease in Lsoftx because both the luminosity and the stellar winds 
of the most massive stars, increase. After the first 3 Myr of 
evolution, the most massive stars end their lifetimes and ex- 
plode as supernovae. The injection of mechanical energy be- 
gins to be dominated by energy released by supernovae explo- 
sions, as the importance of stellar winds rapidly diminishes. 
During the first 35 Myr, this remains the situation for IB mod- 
els because, for a Salpeter IMF, the supernovae rate declines 
slowly, while there exist stars of sufficient mass (of initial mass 
above 8 Mp) to produce supern ovae ( Cervino & M as-Hesse 
II994t iLeitherer & Heckman|[l995h . In the case of EB models, 
^softx is expected to increase slowly after the first 5 Myr, until 
an equilibrium is reached between the formation and destruction 
of stars that end their lifes as supernovae, i.e., at around 40 Myr 
at solar metallicities, according to th e evolutionary tracks con- 
sidered. |LeiAerer^Heckma^ijT995l) presented the evolution of 
both supernova rate and the injection rate of mechanical energy 
for longer term evolution, up to ages of 300 Myr Their Fig. 56 
showed that the asymptotic rate of energy injection is within 0.05 
dex of the value predicted at 30 Myr. Although the computation 
of the soft X-ray luminosities in the CMHK02 models is sim- 
pl istic, th e predictions are in good agreement with the results 
of l Strickl and & Stevens ( 1999), which were computed using hy- 
drodynamical simulations of a young superbubble driven by a 
cluster of massive stars. 

Figure |2] shows the predictions for the Lsoftx /^^fir ratio evo- 
lution with time. It can be seen that Lsoftx/i^FiR increases contin- 
uously with time after an instantaneous burst, even after the first 
5 Myr, due to the fact that Lsottx remains almost constant while 
LpiR decreases rapidly. In extended burst models, the rate of in- 



log [L^/LpJ 

Fig. 4. Lsoftx /i'FiR histograms for the samples of star- 
forming gala xies from Ranalli, Comastri. & Sett 3 (I2003h 
(dashed line), Tiillmann et al. (2006b) (thick solid line) and 
iRosa-Gonzalez et al. ( 2007) (thin solid line). The bins have 
been computed with a slight shift between each other for clarity. 



crease of Lsoftx /Lfir with time, after the first 5 Myr, is smaller. 
As discussed above, we expect the Lsoftx/LpiR ratio to stabilize 
after about 40 Myr in EB models, when an equilibrium between 
formation and destruction of stars susceptible to becoming su- 
pernovae, has been reached. In both cases, there is a rapid in- 
crease in the Lsoftx /Lfir ratio after the onset of the burst (one 
order of magnitude in 3 Myr in IB models), associated with the 
rapid increase in the amount of mechanical energy injected dur- 
ing the first phase of the evolution of the most massive stars. 

In Fig. [3] we plot the evolution of the Lsoftx /Lfir ratio as a 
function of metallicity (Z = Zq and Z - 0.4 xZq) and star forma- 
tion history (instantaneous and extended bursts). Varying metal- 
licity highlights two results: first, low-metallicity stars evolve 
more slowly and have a longer lifetime, such that the evolution 
of LpiR is delayed with respect to solar-metallicity stars. Second, 
the lower the metallicity, the less efficient are the stellar winds, 
and therefore the lower is the amount of mechanical energy re- 
leased into the interstellar medium. As seen in Fig. [3] the net 
effect is a decrease in the value of Lsoftx /Lfir at intermediate 
ages by 0.2 dex, or by up to 0.5 dex within the first 4 Myr of 
evolution. At some ages, the trend is even reversed. For the time 
interval considered, however, the predicted values of Lsoftx /Lfir 
are similar for both values of metallicity. 



3. Observational data sample 

We have compiled Lsoftx and Lfir data for 62 star-forming galax- 
ies, of different types and redshifts, to compare with our model 
predictions. 

The data compilation of iRanalli. Comastri. & Settil (l2003l) 
(hereafte r RC S 3 ) is for star-forming galaxies, from the atlas of 
iHo et al.l (I1997I) . that have detectable X-ray emission in ASCA 
and/or BeppoSAX observations. Only spiral and irregular galax- 
ies from Sa to later types were included in this sample, which 
was complemented by the authors with data of six, well-known 
starburst galaxies observed in t he southern he mis phere. 

The sample compiled by iTiillmann et alj j2006bl) (here- 
after TUL06) is based on th e nine late-type starburst galax- 
ies of iTiillmann et alj (l2006al) (XMM-Newto n data), seven star- 
forming disk galaxies from [Strickland et al.l (|2004b) (Chandra 
data) and seven additional late-type star-forming galaxies taken 
from the literature. The contribution by obvious point-sources 
within the extraction regions was removed from all objects in the 
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Fig. 5. L.softx/i'FiR histogram for the combined sample, over the 
evolution of the Lsoftx/i-FiR computed for IB (dashed lines) and 
EB (solid lines) and fxeft = 0.10 (top) and 0.01 (bottom). The 
horizontal line corresponds to the mean ratio of the complete 
sample, log(Lsofix/LnR) = -3.92. 

original iTiillmann etlTI (l2006ah and [Strickland etal] (l2004bh 
compilations when possible. 

In both samples, the X-ray emission was corrected for 
Galactic neutral Hydrogen absorption, but not for the intrin- 
sic absorption of the galaxies. The published far-infrared fluxes 
were computed u sing the IRAS 60 a nd 100 /urn flux following 
the procedure of Helou et alj (Il988l) . corr esponding to the en- 
ergy emitted within the range 40-120 ^m. lCalzetti et al.l (l2000l) 
found for local starburst galaxies that F(l - I000)/F(40- 120) = 
1 .75 +0.25 . We have therefore multiplied the FIR fluxes provided 
by RCS03 and TUL06 by 1.75, in order to obtain a more realis- 
tic determination of the total amount of energy being reemitted 
in the mid- and far- infrared range. The luminosities were re- 
computed for all objects from the published fluxes by using the 
distances corrected to the Local Group reference frame as given 
in the Nasa Extragalactic Database , assuming Hq - 73 km s"' 
Mpc-i. 

There is some overlap between the two samples. The 
L,oitx/LnR ratio for M82 is -4.08 in RCS03 (based on 
BeppoSAX data) and -3.85 for TUL06 (Chandra). For 
NGC 4631 Lsoftx/ipiR varies between -3.96 (ASCA) and -4.17 
(XMM-Newton), respectively. In these cases we have taken the 
isoftx values provided by TUL06 . 

To compare with star-forming galaxies outside the Local 
Univer se, we consider galaxi es at re dshifts z = 0.01 - 0.67 stud- 
ied by iRosa-Gonzalez et al.1 (l2007h (hereafter ROSA07). Both 
^softx and LpiR data are available for these galaxies, which were 
observed as part of the Chandra Deep Field South (CDFS) 
survey. The X-ray data in thi s sample come from the cat- 
alog of [Alexander et al.l (l2003b . and were not corrected for 
either Galactic or intrinsic absorption by neutral hydrogen. 
Nevertheless, along the line of sight to the CDFS the expected 
Galactic neutral hydrogen absorption in the soft X-rays band is 
only 4.2% (0.02 dex) (Alexander et al. 2003). We exclude ob- 
jects believed to harbour an obscured AGN, or to be dominated 
by low-mass. X-ray binaries (LMXB). In total, we have identi- 
fied 18 objects with detectable soft X-ray, and far-infrared fluxes 
that appear not to be contaminated by either an AGN or LMXB. 

ROSA07 derived the far infrared luminosities for these ob- 
jects in the full 8 - 1000 yum band using Spitz er observations 
at 25 yum. They used the empirical calibration of iTakeuchi et al.l 



(l2005h based on the analysis of a large sample of galaxies for 
which fluxes in the four IRAS bands at 12, 25, 60 and 100 
//m were available. The FIR luminosity computed in this way 
should be consistent with the luminosities derived for the RCS03 
and TUL06 samples, extrapolated to the 1 - 1000 fim range. 
ROSA07 assumed Hq = 70 km s"' Mpc"' to convert fluxes into 
luminosities. 

We note that the Lsoftx values are integrated over the 0.5 - 
2.0 keV band by RCS03, 0.2 - 2.0 by ROSA07, and over the 
0.3 - 2.0 keV band by TUL06, while the model predictions were 
computed for the 0.4-2.4 keV band. We have verified that for the 
typical spectral properties of the diffuse gas in these objects, the 
discrepancies associated with the different bandwidths should be 
smaller than 9% (0.04 dex) in any case. 

We have plotted in Fig. |4] the histograms of the Lsoftx /^^fir 
distribution for each sample. It can be seen that while there is 
a significant overlap between them, the star-forming galaxies 
compiled by RCS03 show the smallest dispersion. The TUL06 
galaxies show generally lower Lsoftx /^^fir values than the ob- 
jects compiled by RCS03, while the sample of ROSA07 presents 
the highest Lsoftx/i'FiR ratios. The mean Zo^(Lsoftx/i'FiR) val- 
ues derived for each sample are -3.94 (0.27) (RCS03), -4.34 
(0.48) (TUL06), -3.37 (0.45) (ROSA07) and -3.92 (0.57) for 
the whole compilation. Values within parentheses correspond 
to the cr dispersion of each sample. We have looked for any 
possible correlation between the Lsoftx /^^fir ratios in the TUL06 
galaxies and their morphological type, but there is no clear trend. 
Both Sc4-Sd and SB galaxies in the sample cover a wide range of 
luminosities and Lsottx/^FiR values and furthermore, they cover 
a similar range in luminosities than the galaxies in the RCS03 
sample. There is no obvious reason for the differences between 
the two samples of local star-forming galaxies. 

4. Discussion 

Figure|2]indicates that our models predict a strong dependence of 
the Lsoftx /i'FiR ratio on the star formation history (instantaneous 
or extended), and on the evolutionary state of the star formation 
process. Moreover, the ratio is also strongly dependent on the ef- 
ficiency of the reprocessing of mechanical energy and UV pho- 
tons into soft X-ray and FIR emission, respectively. Therefore, 
we would expect a significant scatter in the Lsoftx /^fir values 
observed in star-forming galaxies. This is indeed what we find 
in Fig.m as discussed above. 

iMas^Hesse & KunthI (Il999h showed that the star formation 
episodes taking place in compact starburst or HII galaxies are 
generally of short duration, and that their properties can be 
better reproduced by evolutionary synthesis models assuming 
(nearly) instantaneous bursts, than by long-lasting, extended- 
in-time, star-formation processes. The generally strong optical 
emission lines in these galaxies constrain the evolutionary state 
of their massive star-formation episodes to ages below 6 Myr, 
typically within 4-5 Myr On the other hand, star formation 
is expected to proceed during long periods of time in the disks 
of late-type spiral galaxies, generally in the form of individual 
bursts at different times. The accumulation of individual star- 
burst episodes along the disks of thes e galaxies mimick a contin- 
uous star formation process (iMas-Hesse & Kunthll 19991) . so that 
extended star formation models might be a better approximation 
to reproduce their spatially-integrated properties. 

In Fig.|5] we have added to the predictions of the Lsoftx /^^fir 
values the histogram corresponding to the whole sample, as 
well as the average observational Lsoftx /i'FiR value. As discussed 
above, there is a degeneracy between exetr and the age of the star 
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Fig. 6. Lsoftx vs. LpiR for the RCS03 (squares), TUL06 (triangles) and ROSA07 (circles) samples. In the left panel we have over- 
plotted the correlation lines corresponding to the Lsoftx/^FiR ratios predicted by IB models for different ages and Cxeff values. The 
right panel shows the corresponding predictions for EB models. 



formation episode, so that it is not possible to discriminate be- 
tween both parameters without additional constraints on the evo- 
lutionary state. Apart from a few irregular galaxies, for which 
an instantaneous model would describe better the star-formation 
episodes they are hosting, most of the galaxies in the sample 
are large spiral galaxies experiencing long-term star formation. 
Figure 5 shows that the mean Lsoftx/i'FiR ratio can be repro- 
duced by relatively young EB models (after about 10 Myr of 
evolution) with moderate efficiencies txeff of about 1%. Galaxies 
with higher Lsoftx/ipiR ratios would correspond to more evolved 
extended starbursts reaching the evolutionary asymptotic phase, 
with fxeif in most cases below 10%. On the other hand, the galax- 
ies with lower Lsoftx/i'FiR values can only be reproduced by very 
young models, after less than 10 Myr of evolution. The galax- 
ies in the sample have been selected in the various compilations 
for being star-forming or starburst galaxies, i.e., galaxies expe- 
riencing a stronger than average episode of star formation. The 
integrated emission in the far infrared and soft X-ray bands of 
some of these galaxies could be dominated by a single but in- 
tense burst of star formation. While star formation proceeds in 
these objects for a long time, these individual, intense bursts are 
not expected to last longer than a few Myr Some of these star- 
bursts could indeed be rather unevolved. We believe that this is 
what we see in Fig.|5] some of the galaxies with low Lsoftx/^^FiR 
ratios could be dominated by a single, intense and relatively un- 
evolved burst of star formation. A deeper study of the individual 
galaxies would be required to confirm this hypothesis, but it is 
beyond the scope of this work. 

iGrimes et al.l (l2005h computed the Lsoftx/^FiR ratios for a 
sample of ultraluminous infrared, starburst and dwarf starburst 
galaxies (the galaxies classified as starburst are already included 
in the TUL06 sample). For their 7 dwarf starburst galaxies (2 are 
already included in the sample by RCS03), they derived a mean 
log(Lso{tx/LpiR) of -4.0, close to the average value found for 
our complete compilation. The Lsoftx/^^FiR ratio of these galax- 
ies (within the range -3.58 to -4.29) is properly reproduced by 
IB models at a pproximately 4-5 Myr, with e^eff within a realistic 
range 1 - 5% dStrickland & Stevenslll999l:ISummers et al.ll200lL 
|2004|) . The mean /og(Lsoftx/i'FiR) value of the 9 ultraluminous 
infrared galaxies (ULIRG) in their sample is -4.5, clearly below 
our average. This indicates that the star formation processes in 



this kind of galaxy might be relatively unevolved and their emis- 
sion dominated by a single, intense episode of star formation. 

In Figure |6l we plot the observational Lsoftx vs. Lfir values 
for the galaxies in the three samples. We include the predictions 
of IB and EB models at different ages, and for fxefr values be- 
tween 1 and 10%. These plots support the main points of our 
previous discussion: the observational Lsoftx vs. LpiR correlation 
in star-forming galaxies can be well-reproduced by evolution- 
ary synthesis models, assuming realistic parameters: age below 
6 Myr for IB cases, and a spread of young to evolved bursts for 
EB models, a high efficiency in the reprocessing of UV photons 
into far infrared emission, and a moderate (1 - 10%) efficiency 
in the heating of the diffuse interstellar gas by the mechanical 
energy released by massive stellar winds. 

We have analyzed the effect of some intrinsic properties of 
the sampled objects, on the dispersion shown in the correla- 
tion plots. First, as noted above, the measured Lsoftx values did 
not i nclude the correctio n for intrinsic neutral hydrogen absorp- 
tion. lKunthetal.1 (119981) measured the column density of neu- 
tral hydrogen in the line of sight to 8 compact starburst galax- 
ies, by fitting their Lyman a profiles, finding values in the range 
loginn) ~ 19-21 cm"^. We have computed that for the typical 
spectral properties of the hot diffuse gas this correction would 
be within 0.4% - 50% (i.e. < 0.18 dex). A second effect would 
be the contamination of the soft X-ray emission by low mass 
X-ray binaries associated to the underlying, older stellar popula- 
tion. ROSA07 estimated that the contribution by LMXB in their 
sample of spiral star-forming galaxies was generally of few per- 
cent, and concluded that the contamination should be negligible 
for galaxies with SFR > 1 Mq ' . In the RCS03 sample only 7 
objects have a value of SFR below 1 MQyr^\ but their average 
/o§(Lsoftx/LFiR) is -3.7 (in any case within -4.2 to -3.2), i.e., 
with no deviation at all from the total average Lsoftx/LnR ratio. 
Therefore, the contamination of Lsoftx by LMXB does not seem 
to be important. Finally, TUL06 removed the contamination by 
point sources before computing the integrated Lsoftx, so that the 
contamination of this sample by LMXB should be negligible. 
An additional effect is related to the relative strength of the star- 
burst emission compared to the galaxy as a whole. While the soft 
X-ray emission would be associated mostly with the starburst re- 
gions, the far infrared luminosity can include a significant con- 
tribution from the rest of the galaxy, so that the Lsoftx /Lfir ra- 
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tios obtained from spatially-integrated measurements would be 
lower than the intrinsic value produced by the starburst itself. A 
detailed morphological analysis of the galaxies in the sample is 
out of the scope of this paper, but we expect that at least in some 
cases the observational Lsoftx/^-FiR value might be contaminated 
by emission not related to the star formation episode. 

Finally, Lfir has been computed assuming almost complete 
reprocessing of UV stellar continuum photons (E(B-V) = 1), 
as discussed in Sect. 2. Lower values of the extinction, of the 
order of E(B-V) - 0.1, would raise the predicted Lsottx/i^FiR 
ratios by up to 0.15 dex. The effective extinction of these ob- 
jects should be in between both extreme values of E(B-V). In 
conclusion, correcting the soft X-ray luminosity from intrinsic 
photoelectric absorption and/or rejecting the far infrared emis- 
sion not associated with the starburst regions themselves, could 
increase the Lsoftx/ipiR ratios for some objects. On the other 
hand, smaller interstellar extinction values would increase the 
predicted Lsoftx/i'FiR ratio by less than 40%. 

We consider if it is possible to derive a calibration that 
would allow the Lsqax luminosity to be used as a tracer of SFR. 
[Strickland & StevensI ([1999) found from detailed hydrodynam- 
ical simulations that superbubbles accelerated by the release of 
mechanical energy in a starburst, would convert on average ap- 
proximat ely 5% of the in put mechanical energy into soft X-ray 
emission. [Summers et alj ( 1200 1ft analyzed Mrk 33, a dwarf star- 
forming galaxy, and concluded that it is dominated by an intense 
burst 5-6 Myr old, and that the rate of injection of mechanical 
energy from the starburst is approximately 1.2x 10"*' erg s"'. The 
soft X-ray luminosity of the central, extended diffuse gas derived 
by these authors is Lsp ftx ~ 2.2 x 10^^ erg s~' , corresponding to 
fxeif ~ 0.018. Similarlv. lSummers et alj (l200 4t) estimated the me- 
chanical energy injection rate from the starburst in NGC 5253 to 
be Lmech = 3.6 X 10"^" erg s"' . The measured thermal X-ray emis- 
sion associated with the starburst was Ljoftx ~ 4 x 10^^ erg s"', 
yielding exeff -0.01. 

We can therefore assume that fxeff is constrained to be in 
the range 1 - 5% for typical star-forming galaxies. Our evo- 
lutionary synthesis models predict for Instantaneous Bursts (of 
age between 3-6 Myr) that Zo§(Lsoftx/i'FiR) is approximately 
equal to log{Lsoftx/LpiR) ~ (-3.1,-4.0), with a central value 
log(Lso{tx/LpiR) ~ -3.5. On the other hand. Extended Burst 
models predict logiL^oftx/LpiR) ~ (-3.0, -3.4) after 30 Myr of 
evolution, when the number of supernova explosions begins to 
stabilize. For less-evolved, extended episodes (of approximately 
10 Myr of continuous star formation) the predicted ratios would 
be lower, within the range log(Lsofix/ Lpir) ~ (-3.5, -4.0). 

These values are close to the average ratio log(Lsofix/ Lpir) 
~ -3.92 derived from our total sample of star-forming galax- 
ies, supporting the use of the soft X-ray emission as a tracer 
of the star formation rate in starburst galaxies. Using the Ljoftx 
values predicted by the models, as shown in Fig.[Tl we can de- 
rive the calibration of the star formation rate (or strength) as 
a function of Lsoftx- For this we have considered the points at 
which logiL^oitx/LpiR) = -3.5 for IB models, log(L^oftx/LnR) 
= -3.1 for EB cases in the asymptotic phase of evolution, and 
logiL^oitx/LpiR) - -3.7 for extended, but non-evolved bursts, 
with ages of approximately 10 Myr, according to the discussion 
above. The calibrations would therefore be: 

SFR (Mq yr-i) = 2 x lO""*' L^oftx (erg s"') (evolved EB) 

SFR (Mo yr-i) = 8 x lO-^i L,of,x (erg s"') (young EB) 

For comparison, the semiemp irical calibration derived by 
iRanalli. Comastri. & Settil (l2003h is 



SFR (Mq yr-i) = 1.1 x lO-^o L,,i,x (erg s"') 

where we have adapted the original calibration to the whole FIR 
emissi on in the 1 - 1000 //m band, to be consistent with lKennicuttI 
(119981) (whose Lpi R cahbration is u sed as a proxy), and have 
scaled the original iKennicutll (11998) mass normalization to the 
range 2 - 12OM0 to be directly comparable with our results. 
These calibrations should be applicable, within the range of va- 
lidity shown in Fig. [T] for galaxies experiencing an extended 
episode of star formation at a constant star formation rate. In 
the case of coeval, instantaneous starbursts, the parameterization 
would be: 

SFS (M0)= 2 x 10-34 ^^^^^^ (gjg g-i) 

where the star formation strength (SFS) is the total mass trans- 
formed into stars at the onset of the burst. 

The semiempirical cahbration proposed by 

iRanalh. Comastri. & Settil (l2003h is therefore in good agree- 
ment with the calibration derived from synthesis models for 
relatively unevolved star formation episodes, but it would tend 
to overestimate the star formation rate for galaxies that have 
been forming massive stars over a long time, for example tens 
or hundreds of millions of years. 

The models do not in principle predict any dependence of 
the Lsoftx/i'FiR ratios on the strength (i.e. total luminosity) of 
the star formation episodes, as proposed by some authors for 
the hard X-ray luminosity. Nevertheless, we want to remark that 
our models do not study the detailed properties of the medium 
surrounding the star-forming regions. The intensity of the star- 
formation episode could, for example, have a direct effect on 
the dust-grain properties of the interstellar medium, exeff may, in 
addition, be a function of the burst intensity. Both effects could 
create a coiTelation between Lsoftx/i-FiR and the star-formation 
burst intensity, although such a correlation is unclear in Fig.|6l 

5. Conclusions 

We have compared the L.softx/i'FiR values measured in a sam- 
ple of 62 star-forming galaxies with the predictions by our evo- 
lutionary synthesis models, aiming to analyze the validity of 
semiempirical and theoretical calibrations of Lsoftx as a star for- 
mation rate estimator. The main results can be summarized as 
follows: 

1. The Lsoftx /i'FiR ratios are strongly dependent on the effi- 
ciency in the conversion of the mechanical energy released 
by the young massive starburst into soft X-ray luminosity, 
by interaction of the stellar winds and supernova ejecta with 
the suiTounding interstellar medium. From theoretical pre- 
dictions and observational data, we expect an exeif value of 
few percent in starburst galaxies. 

2. Lsoftx/^FiR is also dependent on the evolutionary status of the 
star formation episode. It increases rapidly with time during 
the first 5 Myr of evolution of a massive starburst, and shows 
a slower increase afterwards. After a (nearly) instantaneous 
burst of star formation, Lsoftx decreases slower than LpiR, as 
long as there remains a population of massive stars able to 
collapse as supernovae (up to around 35 Myr). 

3. When star formation proceeds at a nearly constant rate dur- 
ing extended periods of time, the Lsoftx /Lfir ratio is expected 
to stabilize after around 40 Myr, when the number of massive 
stars that produce supernovae has reached an equilibrium be- 
tween death and birth of new stars. 

4. The Lsoftx /Lfir values measured for the sample of star- 
forming galaxies are consistent with the predictions by the 
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models under realistic conditions: relatively young and un- 
evolved star formation episodes and txeff values within 1- 
10%. 

5. A calibration of Ljoftx as a star formation rate estima- 
tor, based on the predictions of evolutionary synthesis 
models, has been derived. The calibration proposed by 
iRanalh. Comastri. & Settil (l2003h is consistent with the pre- 
dictions for relatively unevolved, time-extended bursts of 
star formation. 
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